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The many-body dynamics exhibited by living objects include group formation within a population, 
and the non-equilibrium process of attrition between two opposing populations due to competition 
or conflict. We show analytically and numerically that the combination of these two dynamical 
processes generates an attrition duration T whose nonlinear dependence on population asymmetry 
X is in stark contrast to standard mass-action theories. A minority population experiences a longer 
survival time than two equally balanced populations, irrespective of whether the majority population 
adopts such internal grouping or not. Adding a third population with pre-defined group sizes allows 
T{x) to be tailored. Our findings compare favorably to real-world observations. 



Predator-prey systems have been widely studied by 
many disciplines, including physics [1. Outside the few- 
particle limit, mean- field mass action equations such as 
Lotka-Volterra provide a reasonable description of the 
average and steady-state behavior, i.e. dNA{t)/dt = 
f{NA{t),NB{t)) and dNB{t)/dt = g{NA{t), Nsit)) 
where A^a(^) and Nsit) are the A and B popula- 
tion at time t. Such population-level descriptions of 
living systems do not explicitly account for the well- 
known phenomenon of intra-population group (e.g. clus- 
ter) formation [4 , leading to intense debate concern- 
ing the best choice of functional response terms for 
f{NA{t),NB{t)) and g{NA{t), Nsit)) in order to par- 
tially mimic such effects [2 . Analogous mass- action equa- 
tions have been used to model the interesting non- 
equilibrium process of attrition (i.e. reduction in pop- 
ulation size) as a result of competition or conflict be- 
tween two predator populations in colonies of ants, chim- 
panzees, birds, Internet worms, commercial companies 
and humans [3 in the absence of replenishment. The term 
attrition just means that 'beaten' objects become inert 
(i.e. they stop being involved) not that they are necessar- 
ily destroyed. The combined effects of intra-population 
grouping dynamics and inter-population attrition dy- 
namics have received suprisingly little attention^ |5], 
despite the fact that grouping and attrition are so 
widespread [3, 4\ and the fact that their coexisting dy- 
namics generate an intriguing non-equilibrium many 
body problem. 

In this Letter, we consider explicitly the effect of intra- 
population grouping dynamics on the duration T of at- 
trition between two opposing populations A and B. In 
stark contrast to the standard mass-action theories, we 
find that T exhibits a maximum for highly asymmetric 
predator-predator systems in the absence of population 
replenishment (Fig. 1(a)). This non-monotonic depen- 
dence of T on population asymmetry is remarkably insen- 
sitive to whether the majority (i.e. larger) population ex- 
hibits internal grouping or not (Fig. 1(b)). We show how 



T can be manipulated by the addition of a third popula- 
tion which blocks encounters involving smaller fragments. 
Within the physics community, Redner and co-workers 
had considered a related problem of conflict within a 
clustering population in one dimension, and highlighted 
intriguing connections to a more general class of rough- 
ening phenomena in physics [5 . Eguiluz and Zimmerman 
and others had considered an infinite-range one popu- 
lation coalescence-fragmentation model as a simplified 
version of human opinion format ion [6], while Galam and 
others have considered interesting models involving com- 
petition and conflict [7, 8 . Our work focuses instead on 
the consequences of coexisting grouping and attrition on 
the duration T (i.e. the survival time of the minority 
population). 

The intra-population group dynamics in our model are 
driven by essentially the same coalescence-fragmentation 
processes as Ref. [6 , while the inter-population attrition 
process is essentially the same as Ref. [5 (Fig. 1(c)). 
We have checked that our main conclusions are robust 
to a variety of reasonable generalizations (e.g. randomly 
selecting groups independent of group size, or attrition 
beyond a simple cluster subtraction rule [5 , or allowing 
for a limited number of new recruits over time) and to 
a reasonably wide range of parameter space. Our model 
combines the following specific mechanisms: It is well 
documented that groups of objects (e.g. animals, peo- 
ple) may suddenly scatter in all directions (i.e. complete 
fragmentation) when its members sense danger, simply 
out of fear[4 or in order to confuse a predator [4 . (Cu- 
riously, clusters of inanimate objects such as doubly- 
ionized Argon atoms and animal Hox genes, also exhibit 
such complete fragmentation [9 ). Since a sense of danger 
can arise at any time, our model randomly selects a can- 
didate group for fragmentation at each timestep, with 
probability proportional to its size since larger groups 
have more members and hence are increasingly likely to 
spot danger or be themselves spotted [4]. With proba- 
bility ua or ub for groups of type A or B^ the group 



2 



Duration T 



lumerical 
analytic 




■ A: rigid 6: dynamical clustering 

■ A and B: dynamical clustering 
A and S: rigid 



(b) 




Initial /A population N^{^) 



coalescence 




fragnnentation interaction 



FIG. 1: (a) Duration T of attrition (or equivalently, the 
extinction time or survival time of the smaller population) as 
a function of initial A population A^a(O) and fragmentation 
probability va {i^b = 0.3). Solid lines are analytic (Eq. (2)) 
while the surface is numerical simulation. Na{0) + Nb{0) = 
1000. Qualitative features are unchanged by varying ub- (b) 
Black curve: same as Fig. 1(a) with i/^ = = 0.3. Red 
dashed curve: A contains rigid units (size 10) while B features 
internal dynamical grouping (i.e. clustering) as explained in 
text. Green curve: both A and B comprise rigid units (size 
10). (c) Events in our model. Two groups of same type can 
coalesce, e.g. 6 + 4 = 10. Individual groups can fragment, e.g. 
6 ^ 6 X 1. Two groups of opposite type interact, e.g. group 
of size 6 — 4 = 2 remains. 



fragments completely. If it doesn't fragment, a second 
group is randomly selected with probability again pro- 
portional to size, since any subsequent coalescence and 
attrition events will likely be initiated by pairwise in- 
teraction between individual members in the two groups 
and hence the probability will depend on the number of 
members. If the group is of the same type (i.e. A or 
B) the two groups coalesce, mimicking the observation 
that groups may try to build up their size to increase 
their security 'Coalescence' can simply mean that two 
groups act in a coordinated way, not necessarily that they 
are physically joined. If of opposite type, their interac- 
tion leaves an A or 5 group of size \sa — sb\ sa > sb 
or sb > sa respectively (or zero is sa = sb) where sa 
and Sb are A and B group sizes (Fig. 1(c)). Other forms 
of attrition rule (e.g. stochastic) can yield similar re- 
sults. At time populations A and B hence comprise 
nf{t) and nf{t) groups (i.e. clusters) of size 5, where 
T^snf{t) = NA{t) and ^snf{t) = NB{t). Interactions ^ 
are distance-independent as in Ref. [6 since we are in- 
terested in systems where messages can be transmitted 
over arbitrary distances (e.g. modern human communi- 
cations). Bird calls and chimpanzee interactions in com- 
plex tree canope structures can also mimic this setup, as 
may the increasingly longer-range awareness that arises 
in larger animal, fish, bird and insect groups^. 



Figure 1(a) shows our numerical and analytic results 
(Eq. (2)) for the duration T. The initial condition for 
the numerical simulations comprises isolated individuals, 
however the curve is insensitive to these initial conditions 
since the initial group (i.e. cluster) formation times act 
like a small additive term. The present two-population 
numerical implementation is a straightforward general- 
ization of the one-population version discussed in Ref. 
[6 . The excellent agreement suggests that our analytic 
treatment of internal grouping using a time-averaged in- 
teraction term may have wider application within non- 
equilibrium many body systems. Following the model 
mechanisms discussed above, the probability Qab that 
any A cluster is selected and interacts with a B cluster 
is the sum over all s of the probability for an A clus- 
ter of size s to interact with any B cluster, which gives 
(1 - UA)NA{t)NB{t)/[NA{t) + NB{t)]\ The probability 
Qba is similar, with va replaced by ub- After an inter- 
action, A and B are reduced by the size of the smallest 
interacting cluster, whose average value c (i.e. average in- 
teraction size) is well approximated by unity plus a small 
linear correction term 0.2(1 +x)~"'^(l — x){l — i^a){^ — ^b) 
because clusters are generally very small over a large pro- 
portion of T. Note c is not formally the same as the av- 
erage cluster size - in part because interactions do not 
occur at every timestep and the entire system is actually 
time-dependent - but they tend to take on similar val- 
ues. Employing constant c, the populations after i inter- 
actions become Na {t) = Na{0)- ic, Nb {t) = Nb{0)- ic. 
The probability for an interaction between A and B clus- 
ters after i previous interactions is Q{i) = Qab + Qba 
and hence 
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To reduce A^a(^) and NB{t) by c takes l/Q{i) timesteps 
on average. The time to reduce one population to zero 
is the sum of the timesteps required for each interac- 
tion, until the population is eliminated. Supposing B is 
the smaller population, it requires Nb{0)/c interactions 
to eliminate it, hence the final interaction happens after 
Nb{0)/c — 1 previous interactions. The time T to elimi- 
nate the smaller population B is therefore ^ Q~^{i) with 
i running from to ^^^^^ — 1 . Using i = 7+?/;o(n+l), 
where 7 is the Euler-Mascheroni constant and ipo is the 
digamma function, and Xla+i ~ Yli we obtain 

the duration: 
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When A is the smaller population, the form is identical 
but with A and B interchanged. This T expression de- 
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pends only on the initial populations of A and hence 
T = T{x^{u}) for constant where {v} = (i^a^^b)- 
Differentiation yields a maximum T at Xmax — 0.788 for 
= and TV = 10^, independent of {u}. Numerical 
simulations confirm that T{x^{u}) ~ ^i({^})^2(^), with 
ti{{iy}) ^ l/{2 — ua — i^b) exactly as in Eq. (2), thereby 
supporting our use of a constant average interaction size 
c. Two factors therefore determine the duration T: one 
originates from the grouping dynamics within a given 
population (i.e. ^i({^})), while the other originates from 
the asymmetry (i.e. t2{x)). Replacing sums by integrals, 
Eq. (2) can be approximated in the large population 
limit as: 

Na{0) - Nb{0) i^^NBmNAiO) - Nb{0) + c] 
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^Na{0)-Nb{0). 

The peak at x^ax is robust to a variety of model vari- 
ants, and can be understood as follows: When x ~ 0, 
clusters of A and B are abundant and have a reason- 
ably large average size. Interactions between A and B 
clusters are frequent and the attrition per interaction is 
high, hence T is small. As x increases, with A being 
the larger population, an interaction between an A and 
B cluster is increasingly likely to eliminate the B clus- 
ter completely since the A cluster is increasingly likely 
to be the larger cluster. However the interaction rate is 
decreasing rapidly, and T increases overall. For x ^ 1, it 
may take a long time to find a B cluster however there are 
very few to find, hence T becomes smaller. Interestingly, 
the distribution of time-intervals between interactions of 
A and B clusters is approximately exponential for all x, 
except near Xmax where it becomes approximately power- 
law. Note that if the attrition were to end after a given 
fraction of the initial population is eliminated, the same 
qualitative results would still emerge since the theory is 
essentially invariant under overall changes of timescale. 
Figure 1(b) shows the results of A and/or B adopting dif- 
ferent internal grouping. The duration T remains essen- 
tially unchanged if the larger population chooses a static 
internal structure comprising rigid units of a particular 
size. If the smaller population adopts such rigid units, 
T decreases significantly. Hence T is largely dictated by 
the internal group dynamics of the minority population. 
If both A and B are internally rigid, T is small for all x. 

The top red curve in Fig. 2 compares our the- 
ory to empirical results for human conflicts, while the 
lower two blue curves show the mass- action predictions. 
The mass-action equations traditionally used for at- 
trition are^: (1) dNA{t)/dt = -ULNA{t)NB{t) and 
dNB{t)/dt = —ULNA{t)NB{t), called Lanchester's undi- 
rected mass-action model; (2) dNA{t)/dt = —dLNB{t) 
and dNB{t)/dt = —dLNA{t)^ called Lanchester's directed 
mass- action model [3 where ul and dL are constants. 
'Old' wars are blue circles and 'new' wars are red tri- 
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FIG. 2: Duration T of human conflicts, as function of asym- 
metry X between the two opposing military populations. 
x= \Na{0)-Nb{0)\/[Na{0)+Nb{0)]. Data are up to the end 
of 2008, hence final datapoints for the three ongoing wars will 
lie above positions shown, as indicated by arrows. Lower two 
blue lines are the mass-action results. Upper red curve (i.e. 
Eq. (2)) generated using ua = i^b = 0.7 and [Na{0) + Nb{0)] 
fixed (as in (a)). Changing ua and ub changes height of the- 
oretical peak, but leaves qualitative features unchanged. 



angles, with World War II labelled by both since it is a 
natural dividing point. Since TV 1, we take the end- 
point for the undirected mass-action model to correspond 
to reducing the smaller population to one instead of zero, 
thereby avoiding problems with a continuum description 
of NA{t) and NB{t) near zero. Figure 2(b) offers some 
support for a recent hypothesis in the social science do- 
main, distinguishing between old wars in which A and 
B adopt traditional, fairly rigid, military structures, and 
new wars in which B (and possibly A) adopt more fluid 
tactics akin to our model[10 . By contrast, the 'old' wars 
are well described by both the green curve of Fig. 1(b) 
(i.e. rigid armies) and the traditional mass-action theo- 
ries (blue curves), implying that such internal group dy- 
namics were absent in 'old' wars. 

Figure 3(a) shows that the duration T can be manip- 
ulated by adding a third-party population C which can 
block interactions (Fig. 3(b)). For simplicity, we assume 
the Nc members of C are permanently arranged into nc 
groups each with sc permanent members. Apart from 
peacekeepers in human conflict, C could mimic the tar- 
geted blocking of interactions between particular physi- 
ological clusters. A and B undergo dynamical clustering 
as before, except that if a C group is selected and it 
is bigger or equal to the size of the A and B clusters, 
the interaction is blocked and the two A and B clusters 
are permanently pacified (i.e. neutralized). Figure 2(a) 
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FIG. 3: (a) Black curve as in Figs. 1(a), 1(b), 2(b) with A and 
B undergoing internal dynamical clustering, i/^ = z^b = 0.3. 
Red dashed curve: nc = 100 third-party groups, each of size 
sc = 1- Green dotted curve: nc = 1 third-party group, of 
size sc — 100. (b) Third-party blocking event. If neither A 
or B clusters are bigger than C cluster, then C cluster blocks 
the interaction and permanently neutralizes both clusters. 



shows that if C comprises only a few, large groups (e.g. 
green dotted curve) then T decreases irrespective of the 
asymmetry. Having a few, large C groups means that 
some sizeable battles can be blocked, however it also al- 
lows the build-up of sizeable groups of both A and B 
which in turn makes the typical size of interactions big- 
ger. By contrast, if C comprises many small groups (e.g. 
red dashed curve) T can be much larger, showing a huge 
increase around x ^ 0. If real-time management of the 
C population is possible, this duration profile can be ma- 
nipulated even further. 

We stress that our model findings are not a simple 
consequence of either dilution leading to reaction slow- 
down, or of the specific cluster selection scheme that we 
chose. In our model, as in nature, opposing predator 
groups actively seek each other out at each timestep, 
even if their density is low, making this unlike simple 
chemical dilution, and hence unlike simple mass-action 
equations. Regarding cluster selection, we have verified 
numerically that our main conclusions are unchanged if 
we select clusters independent of size, or use other frag- 
mentation schemes (e.g. binary splitting into two clus- 
ters). This is because the smaller population spends the 
majority of the conflict as very small groups or individ- 
uals, hence the weighting by size is not so important. In 



short, our results emerge from the interplay between pop- 
ulation asymmetry, the presence of clustering, and the 
intentional engagement between the two opposing pop- 
ulations. Although the specific consequences may vary 
by application area, we believe that related phenomena 
lying beyond mass-action predictions will arise in a wide 
range of physical, chemical, biological and social systems, 
whenever intra-population clustering coexists with inter- 
population reactions. 

We thank R. Denney and M. Spagat for very useful 
discussions. 
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